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ABSTRACT 



Context. Some recent observational results impose significant constraints on all the models that have been proposed to explain the 
Galactic black-hole X-ray sources in the hard state. In particular, it has been found that during the hard state of Cyg X-1 the power-law 
photon number spectral index, T, is correlated with the average time lag, < t\-^g >, between hard and soft X-rays. Furthermore, the 
peak frequencies of the four Lorentzians that fit the observed power spectra are correlated with both F and < ?iag >. 
Aims. We have investigated whether our jet model can reproduce these correlations. 

Methods. We performed Monte Carlo simulations of Compton upscattering of soft, accretion-disk photons in the jet and computed 
the time lag between hard and soft photons and the power-law index F of the resulting photon number spectra. 
Results. We demonstrate that our jet model naturally explains the above correlations, with no additional requirements and no addi- 
tional parameters. 

Key words, accretion, accretion disks - black hole physics - radiation mechanisms: non-thermal - methods: statistical - X-rays: stars 



1. Introduction 

The states of Galactic black-hole X-ray sources (GBHs) are 
determined from the following quantities: i) their disk frac- 
tion, which is the ratio of the disk flux to the total flux, 
both unabsorbed, at 2 - 20 keV; ii) the photon number spec- 
tral index, T, of the hard band, power-law component at en- 
ergies below any break or cutoff; iii) the root-mean-square 
(rms) power in the power spectral density (PSD) integrated 
from 0.1-10 Hz; iv) the integrated rms amplitude of any quasi- 
periodic oscillation (QPO) dete cted in the range 0.1-10 Hz 
ilRemillard & McClintockl |2006|) . In particular, the so-called 
hard state is characterized by a power-law X-ray spectrum with 
photon number spectral index 1.5 <r < 2.1, which suffers an ex- 
ponential cutoff at ~ 100 or few 100 keV, and contributes more 
than 80% of the 2-20 keV flux. 

Several models have succeeded i n reproducing the observed 



X-ray sp ectrum in this state (see e.g. Titarchuklll994 ; Esin et al. 
,1997; Po utanen & Fabianlll999l;lReig et al.ll2003l;lGiannios et"aL 
2004). When a GBH is in the hard state, radio emission is 
also detected and a jet is either seen or inferred (Fender 20011). 
Again, several models are successful in reproducing the en- 
ergy spectrum from the radio domain to the hard X-rays (see 
e.g.Markoffet al. 2001; Vadawale et ani200ll:ICorbel & Fended 
12002 ; Marko ff^et al.1 120031; IGianniosI 12005). The multiplicity of 
models that can fit the time average spectrum of GBHs well in- 
dicates that this alone is not enough to distinguish the most real- 
istic among them. 

Timing data of the observed intensity in various energy 
bands provide a totally different "dimension" than the spectral 
one, and can in principle provide information on the location. 
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size, physical conditions and kinematics of the plasma responsi- 
ble for the X-ray emission in these objects. For example, the X- 
ray light curves of GBHs in different energy bands show similar 
variations, but those in the higher energy band generally lag the 
variations detected in the lower energy band JNowak et al.lll999t 
iFordet anil999.) . The time delay, tug, usually depends on fre- 
quency, V. For Cyg X- 1 the time delays decrease with increasing 
frequency roughly as flag oc v'^, where yS =s 0.7. As constrain- 
ing as it is, this re lation has also been expla i ned in more than 
one way (see e.g. Poutanen & FabianI 119991: iReig et all 120031 ; 
iKording & Falckell2004 . 

Even more constraining to models is the fact that the high- 
frequen cy power spectrum flattens as the photon energy in- 
creases (iNowak et alj|1999l) . Equivalent to this is the fact that 
the width of the temporal autocorrelation function of the light 
curves of Cyg X -1 decreases with increasing photon energy 
(iMaccarone et al.l [2000). These observational facts are con- 
trary to what one expects from simple Comptonization mod- 
els, yet possible explanations for them have also been offered 
(iKotov et al.ll200ll;lGiannios et al.ll2004 . 

Significant advances in our understanding of how the X-ray 
mechanism in GBHs (and compact accreting objects in general) 
works can be achieved when we combine variability and spec- 
tral information. This has been possible in the last few years, 
due mainly to R XTE. For example, in the case of Cyg X-1, 
Pott schmidt et al.l (2003) used 130 RXTE observations between 
1998 and 2001 to study the long-term evolution of the source's 
power and energy spectrum in detail. They used 512 s long 
light curves to estimate the power spectral density (PSD), which 
they fitted with a model that consisted of the sum of multiple 
Lorentzian profiles. They generally achieved good fits using four 
broad Lorentzian components with typical peak frequencies of 
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vi ~ 0.2 Hz, V2 ~ 2 Hz, V3 ~ 6 Hz, and V4 ~ 40 Hz, which vary 
in unison with time: their ratios remain constant despite each 
frequency varies by up to a factor of ~ 5. Up to now, no physi- 
cal explanation of these Lorentzians has been offered. As for the 
energy spectrum, they find that a simple model consisting of a 
power-law component of index F, a multi-temperature disk black 
body and a reflection from neutral material component could fit 
the energy spectrum of the source below 20 keV. 

They present convincing evidence that the timing and spec- 
tral properties of the source are closely linked. In fact, a number 
of the correlations they present impose stringent constraints on 
all the models proposed so far These correlations are: 

a) The observed photon number spectral index, T, correlates 
with the observed average time lag, < fi^g >, between 2-4 keV 
and 8-13 keV, averaged over the 3.2-10 Hz frequency band, 
in the sense that < flag > increases as the spectrum steepens. 
Therefore, a model that can successfully explain the spectral 
slope r should also predict the right amount of time delay be- 
tween these specific energy bands and vice versa: a model that 
can fit well the time lag vs. Fourier frequency relation, should 
also predict the appropriate value of the index F for a given time 
delay. 

b) The observed photon number spectral index, F, 
also correlates with the Lorentzian peak frequencies: the 
frequencies incre ase as the spectru m steepens (see also 
[Shaposhnikov & Titarchugi2006l 5007k _^ 

In summary, the lPottschmidt et al] ( l2003h results show that a 
specific value of F, say F = 2, which is relatively easy to obtain 
with a simple Comptonization model, must also be associated 
with specific values of the Lorentzian peak frequencies and a 
specific value of < flag >! These are challenging relations that 
any model must be able to explain. In this work we investigate 
whether the jet model that we have proposed can meet these con- 
straints under reasonable assumptions. 

We think that the jets in GBHs play a central role in all 
the observed phenomena, not only in the radio emission. We 
have shown that Compton upscattering in the jet of soft pho- 
tons from the accretion disk can explain: 1) the X-ray spectrum 
and t he time-lag dependence on Fourier frequency (Reig et al. 
I2OO3L hereafter Paper I); 2) the narrowing of the autocorrelation 
function and the increase o f the rms amplitude o f variability with 
increasing photon energy dGiannios et alj|2004l hereafter Paper 
II); 3) the energy spectrum from radio to hard X-rays of XTE 
Jll 18-1-480, the only source for which we have s imultaneous ob- 
servations for all energy bands (lGianniosll2005L hereafter Paper 
III). 

In this paper we show that the jet model proposed and de- 
veloped in the above work can explain the F - < flag > relation 
of Pottschmidt et al] (l2003h . with no additional parameters or re- 
quirements, if we let just two model parameters, namely the ra- 
dius at the base of the jet and the optical depth along its axis, vary 
over a reasonable range of values. At the same time, the model 
can also explain the "F - Lorentzian peak frequencies v," rela- 
tion qualitatively if we identify l/v,, / = 1,2,3,4, with quanti- 
ties proportional to the characteristic timescale MouJMout, where 
Mout is the rate of mass outflow in the jet and Mout the available 
mass to power the jet. Finally, the model parameter variations 
that are needed to explain the lPottschmidt et al.l (l2003l) relations 
also predict that, on long time scales, the radio variations should 
also correlate with F, which, as we show, is indeed the case in 
CygX-1. 

In Sect. 2 we describe briefly our model, in Sect. 3 we 
present and discuss our results, and we close with our conclu- 
sions in Sect. 4. 



2. The model 

The jet model used here is the same as the one described in 
Papers I, II and III. A summary is given below so that the reader 
becomes familiar with the parameters involved. First, we dis- 
cuss the asumptions made in the model; then we describe how 
our Monte Carlo code works and finally we explain the meaning 
of the various parameters involved. 



2. 1 . Description of tine jet 

We assume that the jet is accelerated close to its launching area 
and maintains a constant (mildly relativistic) bulk speed, vy, 
thereafter The electron-density profile in the jet is taken to have 
the form 



n^{z) = tiQizolz), 



(1) 



where z is the distance from the black hole perpendicular to the 
accretion disk. The parameters no and zo are the electron density 
and the height of the base of the jet. The Thomson optical depth 
along the axis of the jet is 



T|| = noO-TZQ ln(i//zo), 



(2) 



where H is the extent of the jet along the z axis. 

Let nR^iz) be the cross-sectional area of the jet as a func- 
tion of height z- Then from the continuity equation Mout = 
jiR^{z)mpne{z)v\\, we obtain the radial extent R of the jet as a 
function of height z 

R(z) = Roiz/zo)"^. (3) 

Here Rq = (Mout/-'!'wv"oi^'||)'^^ is the radius of the jet at its base, 
Mout is the mass-outflow rate, nip is the proton mass, and vy is 
the component of the velocity of the electrons parallel to the z 
axis. 

For computational simplicity, we assume that a uniform 
magnetic field exists along the axis of the jet (taken to be the 
z axis). The z-dependence of the magnetic field is dictated by 
magnetic flux conservation along the jet axis to be 



B(z) = Boizo/z), 



(4) 



where Bq is the magnetic field strength at the base of the jet. Also 
for computational simplicity we take the electrons' velocity as 
having two components one parallel (vy) and one perpendicular 
(vj^) to the magnetic field. The magnitude and direction of vy are 
fixed; so is the magnitude of Vj^. However, the direction of Vj^ 
is continuously changing, being always tangent to a circle in the 
perpendicular plane with respect to the magnetic field. 

It has been shown (Paper III) that a power-law distribution 
of electron velocities in the rest frame of the flow gives nearly 
identical results. This is because the dominant contribution to 
the scatterings comes from the electrons that have the lowest 
velocity (or the lowest Lorentz y factor), due to the steep power 
law of electron y's required to explain the overall spectrum. 

2.2. The Monte Cario simulations 

For our Monte Carlo simulations we used the same code as in 
Paper I. The code has been checked with an independent code 
that was used in Papers II and III. The procedure we follow is 
described briefly below. 

Photons from the inner part of the accretion disk, in the form 
of a blackbody distribution of characteristic temperature Tbb, are 
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injected at the base of the jet with an upward isotropic distribu- 
tion. Each photon is given a weight equal to unity (equivalently, 
beam of flux unity) when it leaves the accretion disk, and its time 
of flight is set equal to zero. The optical depth T(n) along the 
photon's direction h is computed and a fraction e '^^"' of the pho- 
ton's weight escapes and is recorded. The rest of the weight of 
the photon gets scattered in the jet. If the effective optical depth 
in the jet is significant (i.e., >1), then a progressively smaller and 
smaller weight of the photon experiences more and more scatter- 
ings. When the remaining weight in a photon becomes less than 
a small number (typically 10"^), we start with a new photon. 

The time of flight of a random walking photon (or more ac- 
curately of its remaining weight) gets updated at every scattering 
by adding the last distance traveled divided by the speed of light. 
For the escaping weight along a travel direction, we add an ex- 
tra time of flight outside the jet in order to bring in step all the 
photons (or fractions of them) that escape in a given direction 
from different points on the "surface" of the jet. The longer a 
fractional photon stays in the jet, the more energy it gains on av- 
erage from the circular motion (i.e., Vj.) of the electrons. Such 
Comptonization can occur everywhere in the jet. Yet, a photon 
that random walks high up in the jet has a longer time of flight 
than a photon that random walks near the bottom of the jet. The 
optical depth to electron scattering, the energy shift, and the new 
direction of the photons after scattering are computed using the 
corresponding relativistic expressions. 

If the defining parameters of a photon (position, direction, 
energy, weight, and time of flight) at each stage of its flight are 
computed, then we can determine the spectrum of the radiation 
emerging from the scattering medium and the time of flight of 
each escaping fractional photon. The time of flight of all escap- 
ing fractional photons is recorded in 8192 time bins of 1/128 s 
each. In this way, we can compute the number of photons that are 
emitted by the jet in each time bin and for any energy band. In 
other words, we can create 64 s long light curves for any energy 
band, which correspond to the resulting emission of the jet in re- 
sponse to an instantaneous burst of soft photons that we assume 
enter the jet simultaneously. 

Having created light curves in various energy bands, we can 
now compute delays between any pair of bands. In our case, we 
consider the light curves of the energy bands 2-4 and 8-14 keV 
(identical to the energy ban ds of the light curves that Po ttschmidt 
et al. 2003 used). Following Vaughan & Nowak'('l997^, we com- 
pute the phase lag and, through it, the time lag between the two 
energy bands as a function of Fourier frequency (same as in 
Paper I). Then we compute the average time lag, < flag >, in 
the Fourier frequency range 3.2-10 Hz (like Pottschmidt et al. 
2003). 

2.3. Parameter values 

Our jet model has a number of free parameters, which are the 
same as in Paper I: the temperature Tbb of the soft photon input, 
the extent H of the jet, the radius Rq and the height zq of the base 
of the jet, the Thomson optical depth ry along the axis of the jet, 
the parallel component of the velocity vy of the electrons, and 
the perpendicular velocity component v^. 

We note here that, due to the mildly relativistic flow in the 
jet, the emerging X-ray spectra from the jet depend on the angle 
between the line of sight and the magnetic field direction. In or- 
der to have good statistics in our Monte Carlo results, we com- 
bined all the escaping photons with directional cosine with re- 
spect to the axis of the jet in the range 0.2 < cos 6 < 0.6. Strictly 
speaking, this is inconsistent with our assumption of a uniform 
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Fig. 1. Emergent photon-number spectrum multiplied by E^ 
from our jet model for the reference values of the parameters 
given in Sect. 3. The soft X-ray excess below ~ 2 keV is es- 
sentially the unscattered soft-photon input. When fitted with a 
power law with an exponential cutoff (see Eq. 5), we find F = 2. 1 
and ^cut - 70 keV. 



magnetic field. If the magnetic field and the bulk flow were uni- 
form, then only one value of cos would be appropriate for each 
source. As discussed in Paper II, a range of directional cosines is 
absolutely necessary for reproducing the observed phenomena. 
Therefore, a non-uniform (probably a spiraling) magnetic field 
permeates the jet and this, together with the bulk flow, result in 
the above range of directions. A range of directions is thus im- 
portant, but the actual values of this range are not. 

The values of the model parameters are chosen so that all 
the results of Papers I, II, and III (apart from the radio emission) 
are reproduced, and at the same time the hard band power-law 
component has a slope of F = 2.1 while the average time lag, 
in the range 3.2-10 Hz, between the 2^ keV and 8-14 keV 
bands is 5 ms (as observed by Pottschmidt et al. 2003). We call 
these values the reference values and they are: fcT^bb = 0.2 keV, 
H = lOVy, Ro = 50rg, zo = 5rg, Ty = 3, vy = 0.8c, and v^ = 0.4c, 
where rg = GM/c^ - 1.5 x 10^ cm stands for the gravitational 
radius of a 10 solar-mass black hole. 

In Fig. [1] we show the emerging X-ray spectrum for the ref- 
erence values of the parameters. The spectrum above ~ 2 keV 
can be well-fitted with a power law with an exponential cutoff 



dN/dE oc £-rg-E/E,, 



(5) 



where F = 2. 1 and E^ut - 70 keV. One can also see the soft X-ray 
spectrum (below ~ 2 keV), which is essentially the unscattered 
part of the soft-photon input. 

3. Results and discussion 

3.1. The spectral index - time lag relation 

Our first aim is to investiga te whether we can r eproduce the F 
- < flag > correlation of Pottschmidt et aP (I2003L see their Fig. 
7c) by allowing only a few model parameter values to change by 
a "reasonable amount" around their reference values. 

It was shown in Paper II that keeping the rest of the param- 
eters to their reference values and increasing «o (or equivalently 
the Thomson optical depth along the axis of the jet ry) makes 
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Fig. 2. Correlation between the photon number spectral index F 
and the average time lag < fi^g > between 2-4 keV and 8 - 
14 k eV in the frequency rang e 3.2 - 10 Hz. Empty circles are 
from lPottschmidt et al.l (l2003l) and filled circles are the results of 
our model. The larger grey circle represents the model with the 
reference values of the parameters. 



the emergent spectra harder (the photons are scattered more of- 
ten and gain more energy from the energetic electrons of the 
jet). Similarly, by increasing Rq, the spectrum hardens. This is 
because a larger Rq translates to a larger Thomson optical depth 
perpendicular to the jet axis (we kept the rest of the parameters at 
their reference values) and therefore more photon scatterings. At 
the same time, time lags should also depend on the dimensions 
of the jet and its Thomson depth. 

Thus, since both /?() and no affect the slope of the continuum 
and the time lags, it seems reasonable to run a series of mod- 
els with only /?() and hq varying, while the rest of the parame- 
ters of the model are kept fixed at their reference values. For Rq, 
we considered the values between 12 and 1 10 r^ (with a step of 
AR() - 2), while, for each Ri), ry was left free to vary between 1 
and 10, with a step of Ary = 0.2 (we discuss the model results in 
terms of Ty, rather than no, because it is more informative about 
the average number of scatterings that each photon has suffered). 
For each run, i.e., each pair {Ro, Ty), we derived the photon index 
Tmodei from a fit to the emergent spectrum between 2 and 20 keV 
(i.e., with no exponential cutoff, and thus with a slightly different 
r from that of Eq. 5) and the average time lag, < fiag,modei >, as 
explained in Sect. 2.2. 

In Fig. [21 we repre sent the F- < flag > relation of 
iPottschmidt et al.l (l2003l) (i.e. the observations), while Fig. [3] 
shows the input model parameters {Rq, ry) that generate the 
points (Fmodei, < fiag.modei >) of Fig.|2l We show Only those model 
results that agree best with the observed points. We find that the 
model can reproduce well the observed F- < flag > relation in 
Cyg X-1 if the radius at the base of the jet and the optical depth 
along the jet axis vary between 7?o ~ 5 and 50 Schwarzschild 
radii, r, (rj = 2r^), and ry x; 1.5 - 10, respectively. 

However, even more interesting is that, for the model to ex- 
plain the observational spectral index - time lag relation, Rq and 
Ty must change in a highly correlated way: as Rq increases, ry 
must decrease. In fact, we find that a straight-line fit to the points 
(Fig. [3]) is very close to a t\\ (or no) oc 7?^' relation. Generally 
speaking, when Rq increases, so does the average time lag. The 
spectrum should also harden at the same time. However, if ry (or 
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no) decreases inversely proportional to /?o, it affects the spectrum 
mo re and forces it to s teepen exactly as observed. 

IPottschmidt et al.l (l2003h computed F and < flag > using data 
from RXTE observations that lasted 2 - 7 ks and were performed 
every ~ 10 days over a period of 3.5 years. In effect, these values 
characterize the long-term evolution of the physical characteris- 
tics of the X-ray emitting source. The results that we present in 
this work show that our jet model can explain these long-term 
changes if we assume reasonable variations of just two model 
parameters, namely the width of the jet at its base and the opti- 
cal depth along the jet axis. 

It is important to point out that the no oc 7?y' (or ry oc /J^', 
if all other parameters are kept at their reference values) relation 
implies (via the continuity equation) that Mout(^o) oc /?o, i.e., 
that the mass ejection rate in the jet increases linearly with Rq. 



3.2. The spectral index ■ 
relation 



Lorentzian peak frequencies 



We can also explain, to some extent, the "F-Lorentzian peak fre- 
quencies" relation as follows. As mentioned in the Introduction, 
Mout/Mout is a characteristic timescale for ejecting all the avail- 
able mass. If the ejection process operates as a damped oscillator 
with frequency proportional to the inverse of this time scale, i.e., 
frequency - CMoia/Moia (where C is a proportionality constant), 
then no and Rq, and hence the emitted jet flux, will vary with 
the same frequency. As a result, a Lorentzian with frequency 
CMout/Mout - CnmpV\iR^no/Mout = C'rI/iq will appear in the 
power spectrum. Here we have assumed that Mont and vy are con- 
stant. 

In Fig. |4] we plot F versus v\, i.e. the pea k frequency of 
the first Lorentzian of Potts chmidt et all (|2003) and the F ver- 
sus C'R^no plot that result from our simulations, using the same 
Rq and no parameter values, which reproduce the F— < flag > 
relation of \Pottschmidt et al\ f l2QQ3l) (i.e., the values plotted in 
Fig-Ell for C = 1.2 X 10"^^ cm s~^ Clearly, the model-predicted 
"F-Lorentzian peak frequency vi" relation agrees with the data. 

Here C is one extra model parameter, whose value is diffi- 
cult to interpret at the moment (as it should be associated with 
the unknown details of the ejection mechanism at the base of 
the jet). This parameter is necessary for adjusting the normaliza- 
tion of the F - vi relation. However, according to our model, the 
shape of this relation only relies on the dependence on Rq and no 



Kylafis et al.: Jet model for black-hole sources 



2.6- 



2.4 



2.2 









° C*°s * 



;°CA°oO 



0.1 



Vj (Hz) 



Fig. 4. Correlation between the photon number spectral index 
r and the frequency vi of the first Lorentzi an component of 
the po wer spectrum. Empty circles are from iPottschmidt et alj 
( l2003h and filled circles are the results of our model. The larger 
grey circle represents the model with the reference values of the 
parameters. 



through the combination /?q«(). We believe that the similarity be- 
tween the shape of the observed and the model-predicted F - vi 
relation, when Rq and no take the same values that reproduce the 
F- < flag > relation, is an encouraging result, which supports 
our model. Since the ratios of the Lorentzian peak frequencies 
Vi,i = 1,2,3,4, in Cyg X-1 are constant, there is no need to 
reproduce the F - Vj, j -2,3,4, correlations. 

We note that the v, oc R^rtQ relation, combined with the 
no 0^ ^o' one (which is necessary for the model to explain 
the F - < flag > correlation as we showed above), implies that 
V, oc Rq. In other words, we can explain the relation between 
F and Lorentzian peak frequencies only if these frequencies in- 
crease with increasing Rq. This is contrary to what is usually 
assumed, i.e., that time scales increase with increasing size, and 
hence that frequencies should decrease with increasing Rq. 

In our case, the v, oc Rq relation is due to the fact that we 
assumed that Mout does not depend on radius. If the mass ejec- 
tion rate in the jet increases with Rq (as the model requires so 
that no ^ Rq^, and hence F and < fi^g > correlate as observed) 
then the time scale needed to eject all the available mass to the 
jet should naturally decrease with increasing radius. 



3.3. The F - radio flux relation 

If indeed it is no and Rq that vary with time in such a way that 
no oc /Jq^ then we expect that the radio flux should also vary, 
even if all the other model parameters remain constant. In fact, it 
is reasonable to believe that there will be a correlation between 
this flux and F, time lags, and Lorentzian frequencies. It is this 
possibility that we investigate in this section. 

The radio emission of the jet in our model is a result of the 
synchrotron emission of the fast-moving electrons in the mag- 
netic field of the jet. The accurate calculation of the complete 
synchrotron spectrum can be achieved numerically by dividing 
the jet in slices across the z-axis, and calculating the synchrotron 
emission and absorption in each slice assuming a power-law 
electron distribution (Paper III). There it is shown that the radio 
emission from the jet is characterized by a flat spectrum with 
flux density close to observations. For the purposes of this work. 



we derive the scaling of the radio flux density with the model 
parameters. 

The observed flux density Py that comes for height z of the 
jet scales as (e.g. Rybicki & Lightman 1979) 



Py^C„B'^j^^^yv^RHz)z, 



(6) 



where y„„„ stands for the minimum Lorentz factor of the electron 
distribution, and a for its power-law exponent (Paper III). The 
function €„ depends only on the index a. Using Eqs. (1), (3), 
and (4), the last expression can be rewritten as 



Py <x cX' r^>ov'?.R^Zo(z/zo)' 



(7) 



The synchrotron radiation is strongly self-absorbed below a 
characteristic frequency (the turn-over frequency) Vf. The syn- 
chrotron absorption coefficient is given by the expression (e.g. 
Rybicki & Lightman 1979) 



Aan^rjB' 



(8) 



where A„ only depends on the index a. The turn-over frequency 
in a slice of the jet can be estimated as the frequency for which 
the optical depth to synchrotron absorption becomes unity across 
the width of the jet, i.e., ayR{z) - 1. Solving this expression for 
Zy/zo and using Eqs. (1), (3), (4), and (8), we have 



ZvM)=A>-Sin-r-ViJo""^ 



(9) 



Most of the emission seen at frequency v comes from the 
height of the jet that corresponds to the thick-thin transition Zy. 
Inserting expression (9) into expression (7) and using Eq. (2), we 
have for the radio flux density, Py, that one measures at a certain 
frequency, v. 



fy oc u„K„ Zq Bq T|| y^^^ V 



(10) 



where C'^ is a decreasing function of a. This is expected since 
a steeper electron distribution (with the rest of the parameters 
fixed) leads to weaker radio emission. 

If only Ro and ry vary with time according to the relation 
T|| oc /Jy' as we argued above, then the above equation implies 

n+7 -4 

that Py oc /Jy+'TJg*' or Py cc Ro. In other words, we expect the 
radio flux to increase with increasing width of the jet. Since as 
Ro increases (and no oc R^^) the spectrum steepens, we should 
also expect a positive correlation between Py and F: the radio 
flux should increase as the spectrum steepens. 

In Fig.|5]we plot the average 15 GHz radio flux, Pisghz, as 
observed with the Ryle telescope, binned over a period of +2 
days before and a fter th e time of the RXTE observations that 
IPottschmidt et al.l (l2003h used to determine F. (The results do 
not change significantly if we bin the radio light curve using a 
smaller or a larger bin size.) One can see that there is a general 
trend for i^iscHz to increase with increasing F, as expected, until 
F reaches the value of ~ 2.2. When we use the (Pisghz, T < 2.2) 
data points to compute Kendall's t, we find that t - 0.21. The 
probability of obtaining this value in the case of no association 
between Pisghz and F is less than 0.2%. We conclude that, as 
long as F < 2.2, there is a weak, but statistically significant, 
correlation between i^iscHz and F in the sense that, on average, 
the radio flux increases as the spectrum steepens. 

According to Eq. (10), the break of the Pisghz - T correla- 
tion at F X! 2.2, as the source moves into intermediate/transition 
states (and then to the thermal state), suggests that, as Ro keeps 
increasing, then Bo may decrease. Alternatively, the transition 
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Fig. 5. The 15 GHz flux vs the photon number spectral index, T, 
of Cyg X-1 between 1998 and 2001. 

may be characterized by weaker particle acceleration leading to 
either lower values of ymin or higher values of the electron index 
a. If real, this result offers a hint as to what happens as the source 
approaches the thermal state: the width of the jet increases and 
the magnetic field diminishes (which explains the silence in ra- 
dio) and there is weaker particle acceleration taking place in the 
jet (which explains the weaker X-ray power-law emission). 



4. Summary and conclusions 

In a series of three papers (Papers I, II, and III) it has been shown 
that our jet model can explain a number of observations regard- 
ing black-hole X-ray binaries in the hard state. These are: a) the 
energy spectrum from radio to hard X-rays; b) the time-lag of the 
hard photons relative to the softer ones as a function of Fourier 
frequency; c) the flattening of the power spectra at high frequen- 
cies with increasing photon energy and the narrowing of the au- 
tocorrelation function. 

In this paper we investigate whether our model can explain 
the long term variations in the spectral and timing properties 
of Cyg X-1, which are closely linked. The main result of our 
work is that we can rep roduce the T- < fi^g > correlation of 
iPottschmidt et al.l (l2003l) if we assume reasonable variations in 
just two model parameters around their mean values: the width 
of the jet at its base and the optical depth along the jet's axis 
(or equivalently the electron density at the base of the jet). We 
find that, if ry s; 1.5 - 10, T^o ~ 5 - 50r5 and Ty oc R-^\ then 
the model pre dicts a F- < flag > relatio n that is almost identical 
to the one that 'Pottsch midt et al.l (l2003h observed. Furthermore, 
if the flux variability operates on a time scale which is pro- 
portional to Mout/A^out (i-e. the characteristic time needed to 
eject all the available matter to the jet), then the same t and 
Ro variations can also explain the shape of t he "spectral index 
- Lorentzian peak frequencies" correlation of Potts chmidt et all 
(1200 3). Finally, given the model parameter variations that are 
needed to explain the observed "spectral - timing" correlations, 
the model predicts a "radio flux - F" correlation, which, as we 
show, does indeed hold in Cyg X-1. 

The T|| cc /Jy' relation, which is necessary to explain the ob- 
served F- < flag > relation, is a natural outcome of the model 
if simply Mout '^ Ro- We propose that this is the underlying 
"fundamental" relation that governs the long-term evolution of 
the physical characteristics of the jet, hence of the observational 
spectral and timing properties of the source. The physical reason 
for this relation is unclear at the moment; but whatever its ori- 



gin, the relation itself does make sense because as the rate of the 
mass ejected in the jet increases, it seems reasonable to assume 
that the jet will become "bigger", i.e. /?o will also increase. In 
fact, in this way, one can explain reasonably well all the long- 
term observed changes with variations in one fundamental phys- 
ical parameter of the system, namely the accretion rate. As the 
accretion rate increases, Mout should also increase. If the jet size 
increases proportionally to Mout, then no "^ ^o', ^nd hence the F 
- < flag > relation. We hope that this "accretion rate - ejection 
rate -jet size" relation will provide an important and useful hint 
for all models that try to explain the connection between accre- 
tion flows and outflows, hence the formation of jets in regions 
close to the central engines in accreting compact objects. 

We find that the Mout °c ^o relation predicts that the radio flux 
should increase linearly with the photon number spectral index. 
To the best of our knowledge, such a relation has never been 
suggested or even detected either in Cyg X-1 or in any other 
GBH in the past. We show that a linear "radio flux - spectral 
index" relation does exist on long time scales in Cyg X-1. We 
believe that this result strongly supports our view that a) X-rays 
in Cyg X-1 are produced by Compton upscattering in the jet of 
soft photons from the accretion disk and b) the jet's size and opti- 
cal depth evolve with time because Mout °^ Ro- Furthermore, that 
the "radio flux - spectral index" relation breaks at index values 
higher than ~ 2.2, provides interesting hints into what may be 
happening as the source moves to the intermediate and thermal 
state. Most probably the magnetic field decreases and/or there is 
weaker particle acceleration in the jet (i.e. the electron distribu- 
tion index, a, becomes larger) as the accretion rate increases (as 
inferred from the increase in the disk thermal emission). 

The explanation of the short time-scale variability properties 
of Cyg X-1 is significantly more challenging. Our model cannot 
explain why four Lorentzian components exist in the source's 
power spectrum (assuming of course that the Lorentzians are the 
fundamental building blocks of the power spectrum in Cyg X- 1 
and other GBHs). However, it can explain, to some extent, the 
way they vary with the source's photon number spectral index. 
If indeed the long-term evolution of the source is governed by 
variations in Rq and Mout, it is natural to assume that the short- 
term variations are caused by short-term, random variations of 
Mout (and hence of no or Ty, and finally X-ray flux), around its 
mean value. We find that, if these variations operate on a time 
scale that is proportional to the characteristic time needed to 
eject all the available matter to the jet, then the model can explain 
the shape of the "F - Lo rentzian peak frequencies" relations of 
[Pottschmidt et al.l (l2003i) . There is no reason why the model "F- 
inverse time needed to eject the mass to the jet" relation should 
have the same shape as the observed "F- Lorentzian peak fre- 
quencies". Thus we do not find this agreement to be coinciden- 
tal. Instead, it is one more evidence in favor of the jet model and 
the parameters' evolution we present in this work, and we find it 
operates on long time scales in Cyg X-1. 

Our results suggest that the characteristic frequencies seen 
in the power spectrum should increase with increasing size of 
the source, i.e. Rq, a result that goes against current views. 
However, the characteristic frequencies may not scale directly 
with the source size, but with accretion rate, hence Mout. As ar- 
gued above, it seems reasonable to think that, as the accretion 
rate increases, so too the mass ejection rate to the jet. Perhaps, 
then, it is the accretion rate (through Mout) that drives the varia- 
tions in both /?o and on the characteristic time scales. A physical 
explanation as to how this may happen and as to the value (and 
number) of characteristic time scales, requires knowledge of the 
details of the acceleration/ejection process of the jet, and is be- 
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yond the scope of the present work. However, as we mentioned 
above, we hope that the results we present here will help the in- 
vestigation of the proper mechanism responsible for the launch 
of the jets and winds in accreting compact objects. 

We note that our conclusion, presented in Sect. 3.2, that the 
Lorentzian peak frequencies v, increase with the size Rq of the 
base of the jet, holds for a given source (in our case Cyg X-1). 
If one wanted to extend our model to other sources, then the 
mass of the black hole comes in. The argument goes as follows: 
Making the reasonable assumption that Mout '^ Rq and using 
«o oc Rq^, we find that the v, scale as l/Ro- Since Rq is typically 
a few times Kg, we conclude that the Lorentzian peak frequencies 
V; are inversely proportional to the mass of the black hole. 

In closing, we want to point out that the above results and 
conclusions are consistent with the fact that, in the thermal state 
of Cyg X-1, the time lags vers us Fourier frequency are identi- 
cal to those in the hard state dPottschmidt et al.l l2000h . In the 
thermal state, the radio emission is very weak and possibly un- 
detectable, but we think that the jet is there and that it pro- 
duces, by Comptonization, both the time lags versus Fourier fre- 
quency and the steep power-law energy spectrum above ~ 10 
keV. According to our model, the suppression of both the ra- 
dio and the X-ray emission in this state should be caused by a 
magnetic field suppression and a weaker particle acceleration, 
respectively. 
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